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ABSTRACT 

We present initial XMM-Newton EPIC observations of the prototypical starburst 
galaxy M82. The superwind is seen to extend continuously from the starburst region 
to the X-ray emission associated with the Ha 'cap'. We also find evidence for a ridge 
feature, probably associated with the superwind interacting with an Hi streamer. Nar- 
row band images, centred on individual X-ray lines, show differences in morphology, 
with higher energy lines being less spatially extended, and a systematic shift with en- 
ergy in the region of peak emission. Spectral fits with two thermal and one power-law 
component provide a good fit to the spectra of the nuclear and inner wind regions, 
indicating a multiphase superwind. We discuss the implications of these observations 
on our understanding of superwinds. 
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1 INTRODUCTION 

Starburst events in galaxies represent an important phase in 
the evolution of galaxies and ultimately in the generation of 
cosmic structure (Heckman 2002). Starbursts affect galactic 
structure and evolution by heating and metal-enriching the 
ISM via material outflowing from the starburst, and if the 
event is strong enough, removing some of the ISM via a 
galactic superwind. Indeed, if the starburst event is strong 
enough (or the galaxy small enough) a starburst can remove 
the majority of the ISM, with major implications for future 
star- formation (Mac Low & Ferrara 1999). Superwinds can 
have a major impact on their environments, enriching the 
IGM, suppressing galaxy formation and so on. The impact of 
these events may have been particularly important at high 
redshift (Mori, Ferrara & Madau 2002). 

In order to study in detail the impact of superwinds 
on galactic evolution and the IGM it is necessary to look 
at nearby objects, where we have better spatial resolution 
and photon statistics. Nearby superwind galaxies, such as 
NGC253, NGC1569 and M82, are key objects and be- 
cause superwinds are driven by hot shocked gas, originat- 
ing from supernovae or fast stellar winds, then X-rays are 
probably the most important waveband to study super- 
winds. While Chandra is a superb instrument for disentan- 
gling point sources in galaxies, the large collecting area of 
XMM-Newton makes it better for studying extended low 
surface brightness features. Some important Chandra and 



XMM-Newton results on starbursts and superwinds have 
already been presented (eg Pietsch et al. 2001; Strickland 
et al. 2000; Martin, Kobulnicky & Heckman 2002). In this 
paper we present the first XMM-Newton European Photon 
Imaging Camera (EPIC) results for M82. The first M82 re- 
sults from the XMM-Newton Reflection Grating Spectro- 
graph (RGS) have been presented in Read & Stevens (2002). 

M82 is a much studied, nearby starburst and super- 
wind galaxy (D — 3.63 Mpc, Freedman et al. 1994, so that 
1' = 1.06 kpc). M82 is very IR luminous (Rieke et al. 1980), 
contains a substantial population of young supernova rem- 
nants and luminous HII regions (Pedlar et al. 1999), and 
many luminous super star clusters (de Grijs, O'Connell & 
Gallagher 2001), all indicating a strong and ongoing star- 
burst, probably triggered by a close encounter with M81. 

As regards the M82 superwind, extended, extraplanar 
emission is seen at several wavelengths. Ha observations 
show a filamentary structure consistent with a roughly con- 
ical bipolar outflow (Shopbell & Bland-Hawthorn 1998), 
with implied gas velocities along the cone surface of ~ 
600 km s _1 . The implied Ha morphology is roughly cylin- 
drical for a height above the disk of z < 350 pc, flaring 
out at larger z with an opening angle of 30° (McKeith et al. 
1995). An additional structure, termed the Ha 'cap' lies at a 
distance of ~ 11 — 12 kpc to the NW of M82 and is probably 
associated with the superwind (Devine & Bally 1999). 

M82 has been observed by virtually all X-ray satel- 
lites and only a selective overview is given here. EINSTEIN 
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Figure 1. Left: The XMM-Newton view of M82. An adaptively smoothed and exposure-corrected image of the merged data from the 
MOS1, MOS2 and PN instruments in the 0.2 — 10 keV waveband. The extent of the optical galaxy is shown by the ellipse (11.2' X 4.3'). 
Contours are shown to highlight the low surface brightness emission, and increase by factors of 2. Right: The three-colour EPIC image of 
M82. The 0.2 — 0.5 keV , 0.5 — 0.9 keV and 0.9 — 2.0 keV bands are shown by red, green and blue respectively. The main morphological 
features discussed in the paper are indicated. 



first detected extended X-ray emission associated with M82, 
extending for several kpc above the plane and being well 
correlated with the optical filamentary structure (Watson, 
Stanger & Griffiths 1984; Fabbiano 1988). Strickland, Pon- 
man & Stevens (1997), using ROSAT, also studied the ex- 
traplanar emission, extending to z ~ 6 kpc, with asymmet- 
ric morphology and brightness. Fitting spectra to subre- 
gions in the wind, they found that the fitted temperature 
of superwind emission drops from kT ~ 0.6 keV near the 
nucleus to ~ 0.4 keV in the outer wind. Strickland et al. 
(1997) concluded that the superwind X-ray emission was 
most likely coming from shocked clouds within the super- 
wind rather than the superwind itself. Lehnert, Heckman & 
Weaver (1999) discussed an X-ray feature seen by ROSAT 
lying about 11' N of M82, coincident with the Ha 'cap' emis- 
sion, and concluded that this emission was likely produced 
by shock heating, as the superwind encounters a cloud in the 
halo of M82. This highlights a major unresolved issue - the 
origin of the X-ray emission from the superwind - is it due to 
the wind itself or shocked dense clumps within the wind, and 
how much of the superwind material is from the starburst 
and how much is entrained by the outflow? As highlighted 
by Strickland & Stevens (2000) understanding the origin of 
the X-ray emission is important to understanding the mass 
and energy contained in the wind, which in turn is crucial 
to understanding the impact of superwinds on their host 
galaxies and environments. High resolution Chandra obser- 
vations of M82 have provided insight into the behaviour of 
point sources in the central region. The brightest source is 
located away from the galaxy centre and may be an inter- 
mediate mass black-hole or a beamed X-ray binary (Kaaret 
et al. 2001). 



Here we focus on an initial look at the morphology and 
spectra of the superwind. A more detailed analysis will be 
presented in a later paper. 



2 THE XMM-Newton OBSERVATIONS OF M82 

M82 was observed with XMM-Newton on 2001 May 6 (rev- 
olution 258) for 30.9ksec, with the EPIC MOS and PN cam- 
eras in full frame mode with the medium filter. The following 
are preliminary results from the EPIC MOS1, MOS2 and PN 
data. Details of XMM-Newton and the EPIC MOS and PN 
instruments can be found in Jansen et al. (2001), Turner et 
al. (2001) and Striider et al. (2001) respectively. 

The data from the 3 EPIC instruments have been pro- 
cessed with standard procedures in the XMM-Newton SAS 
(Science Analysis System V5.2). Periods of high background 
were filtered out (leaving 30.0ksec of usable data). Only 
events with a pattern number < 12 and flag selection for 
good events were retained, and photons with an energy in 
the range 0.2 — 10 keV . For the PN detector, 'out-of-time' 
events were screened out. The pipeline processing has iden- 
tified over 330 point sources within the EPIC field, though 
many are spurious and associated with the "out of time" 
events. We defer a detailed study of the M82 point source 
population and the associated field to a later paper. 

3 RESULTS: IMAGING THE SUPERWIND 

In Fig. [T] (left panel) we show a combined MOS1, MOS2 
and PN image of M82. The data for each instrument has 
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Figure 2. Left: The M82 X-ray/Ha overlay. The X-ray contours from the combined EPIC MOS and PN image are superimposed on 
a Ha image. The image size is approximately 17' square. Right: The relative intensities (normalised to the peak values) of the X-ray 
(black) and Ha emission (grey) as a function of distance from the M82 nucleus in a slice running northeast (left) to southwest (right) at 
a position angle of 25°. The X-ray source at a[2000.0] = 09 h 55 m 14 s <5[2000.0] = +69°47'32") is visible to the very left of the plot. 



been corrected for periods of high background and out-of- 
time events, mosaiced together and adaptively smoothed. In 
addition to the large number of point sources in the field, 
the main features to note are the bright nuclear region and 
the extended extraplanar emission, predominantly bipolar in 
nature. The superwind emission is asymmetric, extending to 
a greater distance in the N (~ 13') than the S (7'). For the 
northern wind the emission extends nearly to the edge of the 
field of view. The emission is also clearly highly structured. 

There are a few more additional features to note; the en- 
hanced emission in the N, coincident with the Ha 'cap', and 
a clear X-ray bridge connecting this emission to the main su- 
perwind emission. The X-ray emission associated with the 
Ha 'cap' also seems to be structured and contains a point 
source (which could well be a background object). As had 
been suspected, this X-ray emission is now clearly seen to 
be part of the superwind. 

In Fig. □ (right) we show a three colour image of the 
emission from M82. The 3 colour bands (red: 0.2 — 0.5 keV , 
green: 0.5 — 0.9 keV , blue: 0.9 — 2.0 keV ) have been chosen 
to highlight the softer superwind emission. In addition to 
some harder point sources in the field, much of the super- 
wind emission is of a fairly consistent colour, implying little 
temperature change outside the inner region. There does ap- 
pear to be a softening towards the edge of the superwind. 
As can be seen from the spectral fits discussed later we see 
a reduction in both the fitted temperatures and in the level 
of absorption suffered. Both these effects are responsible for 
the softening of the wind emission. 

The X-ray emission associated with the Ha 'cap' is con- 
centrated in the 0.5 — 0.9 keV (green) waveband. Another 
interesting morphological feature is a ridge of emission in the 
SE portion of the wind (and which is indicated in Fig. 0. 
Similar ridge like features in the X-ray emission from super- 
winds have been seen in NGC 253 for example (Strickland 
et al. 2000). In the N wind there is a region with a deficit 
of emission, which we term the 'dark lane', which is well 
correlated with HI emission (see later). 



In Fig. [5] we show X-ray contours from the combined 
EPIC MOS1, MOS2 and PN image, superimposed on a Ha 
image of M82. Note that this Ha image does not extend out 
to the Ha 'cap'. There is a good general tie-up in the super- 
wind, between the filamentary Ha and the X-ray emission. 
In the northern wind there are indications of the X-ray emis- 
sion more closely following the filamentary structure, with 
Ha emission associated with the tongue of X-ray material 
connecting to the Ha 'cap' emission. 

Fig- |U also shows how the relative intensities of the X- 
ray and the Ha emission vary with distance, in a NW-to- 
SE slice (at a position angle of 25°) through the nucleus of 
the galaxy. Both profiles have been background-subtracted. 
Though the Ha emission appears more centrally concen- 
trated than the X-ray, several similarities in the profiles are 
evident. Both profiles show, that within the inner wind (up 
to 100") the southern side is brighter. This is likely due to 
the fact that, due to the orientation of the system, we are 
able to see the southern side of the inner wind far more 
clearly than the obscured northern side. Beyond this how- 
ever, the situation reverses; both the X-ray and the Ha pro- 
files show a flattening of the intensity distribution to the N, 
whereas to the S, both profiles fall more rapidly, and disap- 
pear into the background. The northern wind is larger and 
extends further than the southern wind. 

There are some striking correlations (and anti- 
correlations) between the X-ray emission and HI features. 
Contours of HI emission, constructed using the HI data 
of Yun, Ho & Lo (1994), are shown superimposed on the 
XMM-Newton image of M82 in Fig. El An apparent 'hole' 
in the HI emission is visible to the NW, coincident with the 
northern X-ray wind. It is likely that the hot X-ray wind, 
in travelling out into the IGM, has blown this hole in the 
intervening cold neutral material. Note though, that the HI 
features have masses comparable to those implied for the 
superwind, and so could provide a substantial obstacle for 
the outflow (Yun et al. 1993). The large HI feature to the 
east of the northern wind appears to be sharply bordered 
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Figure 3. Contours of HI emission (constructed using the HI 
data of Yun et al. 1994) superimposed on the combined EPIC 
XMM-Newton image of M82. 




Figure 4. Narrow band images, centred on the following lines: 
OVIII (top left), FeXVII (0.82 keV ) (top right), NeX (bottom 
left) and SiXIII (bottom right). Each image is 12' to a side, and 
the contours are at the same levels in each image and increase by 
factors of two. The axes are located at the centroid of the SiXIII 
emission and show the systematic shift in peak emission. 



by X-ray emission from the superwind, and may be colli- 
mating the flow. The 'dark lane' in the X-ray emission in 
the N wind seems also to be associated with enhanced HI 
emission. The 'dark lane' structure seen in the X-ray image 
could be due to either absorption by foreground HI emission 
or due to an interaction is not clear. The fact that we do not 
see a hardening of the spectra in this region in the 3 colour 
diagram (Fig. suggests that it is not simply foreground 
absorption. 

The southern wind does not show such an obvious HI 
'hole' as the northern wind, but the X-ray ridge is inter- 
esting, lying at the southern edge of the large HI streamer 
extending to the east. The HI streamers to the southwest 
(in the direction of M81) do not show any significantly en- 
hanced X-ray emission. Also, the Ha 'cap' does not show 
any associated HI emission, and hence has been assumed to 
be due to an ionized cloud. 

The large effective area of the EPIC instruments on 
XMM-Newton allows the generation of narrow band images, 
centred on particular lines, which gives further information 
on the thermal and (potentially) the abundance structure 
of superwinds. In Fig. 0]we show 4 such images, with each 
image centred in energy space on 4 lines, from low energy 
lines (OVIII, 0.65 keV) up to higher energy lines (SiXIII, 
1.85 keV). The differences in morphology are marked, with 
the spatial extent of the higher energy lines being much 
smaller than the lower energy lines. There is also a system- 
atic shift in the peak of the emission, with the emission peak 
of the lower energy lines being shifted to the S and SE by 
up to ~ 0.5'. The orientation of the M82 disk is such that 
the NW edge of the disk lies nearest to us (McKeith et al. 
1995). As such, we are able to see the central regions of the 
X-ray outflow much more easily in the S than in the N. The 
energy shift in peak emission is very likely due to a combi- 
nation of two effects — the central plume of the very inner 
starburst wind (visible only in the south) above (i.e. to the 



SE of) the nucleus is visible in the lower energy band images, 
whereas the actual starburst nucleus is completely absorbed 
by the intervening disk of M82. As we move towards higher 
energies, the contribution of the central plume of the very 
inner wind emission decreases, whilst the starburst nucleus 
becomes progressively less absorbed, and more visible. An 
almost identical situation is seen in the other famous nearby 
starburst galaxy NGC 253 (Pietsch et al. 2001). Note also, 
that the X-ray spatial structure in the central regions of 
M82 has also been studied by analysing the dispersion and 
cross-dispersion profiles of the X-ray spectral lines observed 
with the RGS (Read & Stevens 2002). These results agree 
very well with the results presented here. 



4 RESULTS: X-RAY SPECTRA 

Previous X-ray spectral analyses of M82 have shown that 
multiple spectral components are needed (for instance, 
Moran & Lehnert 1997; Ptak et al. 1997), indicating the 
presence of multi-phase X-ray emitting material. In the light 
of these results we adopt a spectral model with two ther- 
mal and one power-law components (though we do note the 
ambiguities in results and physical interpretation that can 
result from such models - see for example, Dahlem et al. 
2000). In addition, we find that one and two component 
models give substantially inferior fits. 

In this first EPIC paper we present a limited number 
of results for a few regions of M82 and its superwind. More 
comprehensive results, dealing with the spectral variations 
in the superwind are deferred to a later paper. We have 
extracted spectra from all instruments, defining two con- 
centric source regions, a nuclear region of radius ~ 0.5' and 
an annular inner wind region around the nuclear region, ex- 
tending from an inner radius of 0.5' to an outer radius of 
1.3'. The background has been taken from a circular region 
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Table 1. M82 spectral fits for the nuclear region, inner wind region and the Ho 'cap' region. The 90% confidence intervals are quoted, as is 
the emission fraction from each spectral component. The abundance value applies to both the mekal components for the nuclear and inner 
wind region. The model component normalisations are also given — for the mekal models the normalisation is (10~ 14 /47r.D 2 ) j n^n^dV , 
with V the volume and D the distance, while for the power-law models the units for the normalisation are photons keV - 1 cm" 2 s -1 at 
1 keV). The luminosity is the absorbed broad band luminosity of the region, assuming D = 3.63 Mpc. 
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reasonably close to these regions but away from the dif- 
fuse emission. We have used a three component model to fit 
the PN and MOS1 and MOS2 spectra simultaneously, two 
thermal {mekal) components with the abundance parame- 
ter tied together and one power-law component, with each 
of the 3 components having an independent column. The 
physical justification for this model is that in the nuclear 
region, in addition to hot multi-phase superwind gas, we ex- 
pect harder emission from the point source population and 
also possibly inverse Compton emission (Moran & Lehnert 
1997). In the inner wind region, again in addition to the 
hot gas, we expect to detect emission from unresolved point 
sources (background AGN as well as sources associated with 
M82). The best-fit model results, absorbed luminosities etc 
of each component are given in TableQand shown in Fig.|5] 
The approximate number of counts in each of the spectra 
are: nuclear region - PN: 99000, each MOS: 47000; inner 
wind region - PN: 65000, each MOS: 28000 and Ha cap re- 
gion - PN:4650, each MOST800. The power-law component 
in both regions is highly absorbed (though the values are 
not inconsistent with the HI values for absorption in front 
of SNRs in M82 quoted by Wills, Pedlar & Muxlow 1998) 
and has a slope broadly comparable to that of X-ray bina- 
ries and AGN. Both thermal components are relatively soft, 
with values broadly comparable to those seen with previous 
instruments. The single metallicity value is (apparently) rea- 
sonably well constrained at ~ 0.3^0, though the reliability 
of this must be questioned (see below). The X-ray emission 
spectrum of the Ha 'cap' region is also shown in Fig. |K| On 
account of the relatively low count statistics the spectra are 
fitted with a single temperature mekal model (see Table 0, 
with kT ~ 0.6 keV . No firm constraints can be placed on 
the metallicity of this region. 

From these fits we can estimate masses of the X-ray 
emitting gas in these regions (assuming spherical symme- 
try). Summing the masses in the two thermal components, 
we find wind masses of 1.5 x 10 7 r/ 1 ^ 2 Mq for the nuclear re- 
gion and 1.7 x 10 1 n 1 ' 2 Mq for the inner wind region, where 
r\ is the volume filling factor of the X-ray emitting material 
(see Strickland & Stevens 2000 for a discussion of rj). These 
values are broadly comparable with those of Strickland et 
al. (1997), given the different assumed geometries. 

We expect the thermal emission from the superwind to 



come from both nuclear processed material expelled from the 
starburst and material entrained into the superwind from 
the galactic halo. The measured wind abundances can then, 
in principle, be an important constraint on superwind ejec- 
tion processes, and be used to estimate the level of heavy 
element pollution of the IGM by the starburst. We can in- 
vestigate the abundances of the nuclear and inner wind re- 
gion by grouping elements into 3 groups: (1) He, C, and N 
(abundance set to solar), (2) Fe group - Fe and Ni (asso- 
ciated with SNe I), and (3) a elements - O, Ne, Mg, Si, 
S, Ar and Ca (associated with SNell). Fitting again with a 
two temperature plus power-law model, we find that (with 
the abundances tied between the two models) that the Fe 
group abundance is fitted to be O.43o;to^(Fe,0) and for the 
a group O.37o;35.Z( a! 0), implying a slightly subsolar a/Fe ra- 
tio. For the inner wind region the corresponding values are 
O.29o;|ti^(Fe,0) and 0.54o;|oZ( aiQ ) , implying a supersolar ra- 
tio (cf. Martin et al. 2002 results for NGC 1569). 

The presence of an Fe line at ~ 6.5 keV in the M82 
spectrum has been suggested by BeppoSAX observations 
(Cappi et al. 1999). To verify this we have fitted a power- 
law+Gaussian model to the 5—10 keV spectra of the nuclear 
and inner wind regions. For the nuclear region, we find a 
line at 6.6l|:|| keV, width a = O.lSgloJ keV and a line flux 
of 2.2f ; J x 10 -5 ph cm -2 s _1 . For the inner wind region, 
we find a line at 6.43^|3 keV and a line flux of 6.83 8 q 9 x 
10" 6 ph cm s , with a poorly constrained. These values 
are broadly consistent with Cappi et al. (1999). 



5 SUMMARY AND CONCLUSIONS 

In this paper we have presented initial results from 
XMM-Newton EPIC observations of M82, which show the 
extraplanar X-ray emission associated with the superwind 
of M82 in greater clarity than has been seen before. The 
main features are that the superwind emission extends out 
to a height of at least 14 kpc in the N and 7.5 kpc in the S 
and is highly structured. The superwind is continuous out to 
the X-ray emission associated with the Ha 'cap'. The X-ray 
and Ha profiles are similar, with the S wind brighter but 
overall less extended. The are several interesting X-ray /HI 
correlations. In the N wind the superwind appears to have 
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blown a hole in the HI gas, as revealed by a depression in 
the HI emission coincident with a region of X-ray emission. 
This penetration may explain why the N wind is more ex- 
tended. The is also some evidence of the superwind being 
collimated by HI gas. An X-ray ridge-like feature towards 
the SE edge of the wind may be due to an interaction with 
an HI streamer. 

The X-ray spectra from the central regions of M82 can 
be fitted with a three component model (2 thermal+power- 
law) with thermal components with kT = 0.5 keV and 
0.9 keV (nuclear region), and kT = 0.4 keV and 0.7 keV 
(inner wind region) . The emission is consistent with a model 
with multiphase thermal emission from the superwind, along 
with power-low emission from compact sources either in the 
galaxy or background AGN. 

In summary, we have presented only a short overview of 
the XMM-Newton EPIC data on M82. More thorough work, 
which would include dealing with the point source contam- 
ination of the superwind, is needed for a detailed analysis 
of the superwind. The X-ray emission from the superwind 
is extremely complex and these observations have revealed 
new evidence of interactions of the superwind with circum- 
galactic material. Future work will concentrate on different 
aspects, such as the bright central point sources, the spatial 
and thermal structure of the superwind and the relationship 
to the observed Ha emission and the overall energetics of the 
starburst and the superwind. 
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